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DESIGN  OF  INTERCHANGER  FOR  DETERMINING  THE  RATIO  OF  C  MIXED  TO  C 

P  P 

UNMIXED  IN  THE  HELIUM -NITROGEN  SYSTEM 

By 

Robert  E.  Barieau— ^ 

ABSTRACT 

In  this  report  the  design  parameters  for  an  interchanger  to  be 

used  in  determining  the  ratio  of  C  mixed  to  C  unmixed  in  a  helium- 

P  P 

nitrogen  system  have  been  calculated.  It  was  found  that  an  inter¬ 
changer  constructed  of  3/16  inch  O.D.  hard -drawn  copper  tubing  would  be 
satisfactory  for  a  design  pressure  of  5000  psia.  An  arrangement  of 
5  helium  tubes,  19  nitrogen  tubes,  and  31  mixture  tubes  will  provide 
the  required  flow  rate  with  sufficiently  low  pressure  drops  through 
the  interchanger  (maximum  Ap  =  1.7  psi) .  The  interchanger  will  pro¬ 
vide  sufficient  heat  transfer  surface  with  a  length  of  10.5  feet. 

Since  it  is  assumed  that  thermocouples  for  temperature  measure¬ 
ment  will  be  placed  on  the  copper  tube-wall  exterior  instead  of  in 
the  flowing  gas  stream,  calculations  were  made  to  determine  the 
additional  interchanger  surface  needed  to  provide  accurate  indication 

1J  Project  Leader,  Thermodynamics,  Helium  Research  Center,  Bureau  of 
Mines,  Amarillo,  Texas 


Work  on  manuscript  completed  October  1964. 
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of  the  final  equilibrium  gas  temperature.  It  was  found  that  less 
than  one  foot  of  additional  length  on  each  end  of  the  interchanger 
would  provide  such  temperature  readings  within  0.01° F. 

Weight  of  solder  required  in  fabrication  has  been  calculated 
and  heat  leak  and  insulation  requirements  have  been  estimated. 

INTRODUCTION 

As  part  of  our  experimental  determinations  of  the  enthalpy  of 
helium-nitrogen  mixtures,  we  will  determine  the  ratio  of  the  heat 
capacity  of  the  mixture  to  that  of  the  same  components  unmixed. 

This  method  has  been  described  in  a  Report  of  Invention  by  Robert  E. 
Barieau,  entitled  "The  Adiabatic  Mixing  Flow  Calorimeter  and  Double 
Heat  Exchanger  Method  for  the  Determination  of  the  Enthalpies  of 
Mixtures",  dated  August  9,  1963.  The  invention  referred  to  consists 
of  the  following:  Unmixed  streams  of  helium  and  nitrogen  are  brought 
to  the  same  temperature  by  means  of  a  high  temperature  heat  sink. 

They  are  then  cooled  by  passing  through  a  heat  exchanger,  pass  to  a 
low- temperature  heat  sink  and  are  cooled  to  its  temperature,  pass 
through  a  second  heat  exchanger,  where  they  are  warmed,  pass  to  the 
high  temperature  heat  sink  where  they  are  brought  to  its  temperature, 
then  pass  to  an  adiabatic  calorimeter  where  the  gases  are  mixed  and 
the  delta  T  of  mixing  is  determined.  The  mixed  stream  then  goes  to 
the  high  temperature  heat  sink,  is  brought  to  its  temperature,  then 
passes  to  the  second  exchanger,  where  it  is  cooled  by  giving  up  heat 
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to  the  unmixed  streams.  The  mixed  stream  then  passes  to  the  low 
temperature  heat  sink  where  it  is  brought  to  its  temperature,  then 
passes  through  the  first  heat  exchanger  where  it  is  warmed  by  the 
unmixed  streams. 

This  report  gives  the  calculated  design  characteristics  of  the 
interchanger s  mentioned  in  the  above  description. 

DESCRIPTION  OF  INTERCHANGER 

This  interchanger  will  be  constructed  of  3/16-inch  O.D.  hard 
drawn  copper  tubing  with  0.042-inch  wall  thickness.  There  will  be 
a  total  of  55  parallel  tubes,  each  10.5  feet  long,  for  a  total  length 
of  577.5  feet.  The  total  weight  of  copper  tubing  will  be  40.85 
pounds.  The  tubes  will  be  soft-soldered  into  a  close-packed  cylindri¬ 
cal  array.  Because  the  tubing  is  to  be  hard  drawn,  and  we  desire  a 
safe  working  pressure  of  5000  psia,  it  will  not  be  possible  to  use 
ordinary  tin- lead  soft  solder  as  this  would  anneal  the  copper  tubing. 
Therefore,  a  bismuth,  lead,  tin,  cadmium,  low  melting  alloy-melting 
below  200° F ,  will  be  used  to  solder  the  copper  tubes  together.  The 
weight  of  low  melting  alloy  will  be  6,8  pounds.  When  packed  into  a 
close-packed  cylindrical  array  the  outside  diameter  of  the  tubing 
will  be  1.46  inches.  The  tubing  will  be  surrounded  by  Linde  CS-5 
insulation,  0.52  inch  thick.  The  outside  diameter  of  the  insulation 
will  be  2.5  inches.  The  weight  of  insulation  will  be  2.6  pounds. 

The  insulation  will  be  surrounded  by  a  vacuum  tight  tube,  as  it  is 
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necessary  to  evacuate  the  insulation  to  about  1  micron  (10  mm)  of 
mercury. 

The  helium  gas  will  be  carried  by  five  of  the  parallel  copper 
tubes.  The  nitrogen  gas  will  be  carried  by  19  of  the  parallel  cop¬ 
per  tubes  and  the  mixture  stream  of  helium  and  nitrogen  will  be 
carried  by  the  remaining  31  parallel  copper  tubes.  The  total  weight 
of  the  interchanger  exclusive  of  the  outside  tube  will  be  50.25 
pounds . 

SUMMARY  OF  CALCULATED  OPERATING  PERFORMANCE  OF  INTERCHANGER 

A  maximum  of  20  standard  cubic  feet  per  minute  of  nitrogen 
will  be  available  from  the  Amarillo  Helium  Plant  at  a  pressure  of 
600  psia.  A  maximum  of  20  standard  cubic  feet  per  minute  of  grade-A 
helium  will  be  available  from  the  Amarillo  Helium  Plant  at  a  pres¬ 
sure  of  2500  psia.  After  mixing,  the  mixture  will  be  returned  to  the 
Amarillo  Helium  Plant,  at  essentially  atmospheric  pressure,  for  recy¬ 
cle  purification. 

We  have  calculated  the  operating  performance  of  this  inter - 
changer  for  mixture  compositions  of  25-75,  50-50,  and  75-25  helium- 
nitrogen  mixtures,  and  also  for  the  pure  components.  The  calcula¬ 
tions  were  made  for  a  pressure  of  500  psia.  The  inlet  temperature 
of  the  pure  helium  and  nitrogen  streams,  at  the  warm  end  of  the 
interchanger,  was  chosen  as  68° F  (20° C)  and  the  inlet  temperature 
of  the  helium-nitrogen  mixture,  at  the  cold  end  of  the  interchanger 
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was  chosen  as  32°  F  (0°  C) . 

The  calculated  pressure  drops  and  temperatures  of  the  streams 
are  given  in  table  1. 

It  is  seen  that  the  AT  at.  the  warm  end  of  the  interchanger  is 

f 

less  than  4° F .  Also,  the  pressure  is  so  low  that  pressure  drop  will 
affect  the  heat  capacity  to  less  than  0.1  percent. 

It  is  seen  from  table  1  that,  the  smallest  rate  of  heat  transfer 
is  for  pure  helium,  amounting  to  506  Btu  hr.  ^ .  Using  Linde  CS-5 
insulation,  I  have  calculated  the  heat  leak  to  be  0.50  Btu  hr.  ^ ,  or 
of  the  order  of  a  0.1  percent  of  the  heat  transferred. 

The  details  of  the  various  calculations  follow. 

SAFE  WORKING  PRESSURE  OF  3/16  INCH  O.D.  COPPER 

TUBING  WITH  0.042  INCH  WALL 

If  t  is  the  wall  thickness  in  inches  and  R  is  the  inside  radius 
in  inches,  then 

t  =  0,042  inch. 

2R  =  3/16  -  2  X  0.042 
2R  =  .1875  -  .084  =  0.1035 
R  =  0.05175  inch 
t/R  =  0.8116. 

Since  t/R  >  1/2,  it  is  necessary,  in  calculating  working  pres¬ 
sures,  to  use  the  formulas  applicable  to  thick  wall  cylindrical  shells. 
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TABLE  1.  -  Calculated 

pressure  drops  and  temperatures  of 

inlet  and 

outlet 

streams 

under  various  relative  flow  rates. 

P  =  500 

ps  ia 

Stream 

Gas 

SCFM 

rjP 

tl 

L  L 

T  -T 

1  2 

Q,Btu  hr  ^ 

Ap , ps i 

1 

He 

20 

35.1065 

68.0000 

3.1065 

-  506 

1.669 

2 

He 

20 

32.0000 

64.8935 

+  506 

0.069* 

3 

* Laminar  Flow 


1 

He 

20 

35.5505 

68.0000 

3.2473 

-  500 

1.669 

2 

He-N 

80/3 

32.0000 

64.7527 

+  748 

0.227 

3 

-2 

20/3 

35.3765 

68.0000 

-  248 

0.098 

1 

He 

20 

35.8399 

68.0000 

3.6602 

-  495 

1.669 

2 

He-N 

40 

32.0000 

64.3398 

+  1224 

0.642 

3 

»2 

20 

36.0805 

68.0000 

-  729 

0,664 

1 

He 

20/3 

35.4078 

68.0000 

3.6570 

-  167 

0.283 

2 

He-N 

80/3 

32.0000 

64.3430 

+  898 

0.397 

3 

»2 

20 

36.0102 

68.0000 

-  731 

0 . 664 

1 

2 

N2 

20 

32.0000 

64.1402 

3.8598 

+  734 

0.282 

3 

N» 

20 

35.8598 

68.0000 

-  734 

0.664 

A  ' 

/ 
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These  formulas  are  (_1) 


(_1)  ASME  Boiler  and  Pressure  Vessel  Code.  Section  VIII.  Unfired 
Pressure  Vessels,  pages  9  and  121. 


(1)  Circumferential  stress 


t 


i);  2 


SE  +  P 
SE  -  P 


(2)  Longitudinal  stress 


i);  z 
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where  t  =  wall  thickness,  inches 

P  =  allowable  working  pressure,  psi 
R  =  inside  radius  of  shell,  inches 
S  =  maximum  allowable  stress  value,  psi 
E  =  joint  efficiency 
(E  =  1  for  tubes) 


with  t/R  =  0.8116, 

Z1/2  =  1.8116 


Z 


3.2819 
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and  for  circumferential  stress 

7  _  1  +  P/S 

1  -  P/S 

Z  -  Z  P/S  =  1  +  P/S 

p/c  =  Z  -  1  _  2.2819 
7  Z  +  1  4.2819 

P/S  =  0.5329 

For  longitudinal  stress,  we  have 

P/S  =  Z  -  1  =  2.2819 

Therefore,  the  determining  expression  is  P/S  =  0.5329.  Copper 
tubing  comes  in  the  following  conditions:  annealed,  light  drawn,  and 
hard  drawn.  Below  100° F ,  the  maximum  allowable  stress  values  for 
these  conditions  are  given  in  table  2. 


TABLE  2.  -  Allowable  stress  and  safe  working  pressure  of  copper  tubing 


Condition 

S 

ps  i 

P  =  .5329  S 
psi 

Annealed 

6  ,000 

3,200 

Light  Drawn 

9,000 

4  ,800 

Hard  Drawn 

11,300 

6,000 

The  maximum  allowable 

working  pressure  is 

also  given  in  table  2. 

For  4500  psi  operation,  the 

copper  tubing  must 

be  either  light  or  hard 

13 


drawn.  For  5000  psi  operation  the  copper  tubing  must  be  hard  drawn. 

We  have  decided  to  specify  that  the  tubing  be  hard  drawn.  We  will 
then  plan  on  a  top  working  pressure  of  4500  psi  and  we  will  set  our 
relief  valves  at  5000  psi.  The  ASME  Boiler  and  Pressure  Vessel 
Code  (2^)  indicates  that  the  maximum  allowable  stress  in  copper  tubing 

(2J)  Reference  1,  page  90. 

\ 

is  the  same  at  400° F  for  annealed,  soft,  and  hard  drawn  copper  tubing. 

This  means  that  heating  hard  drawn  copper  tubing  to  400°  will  anneal 
it.  As  the  melting  point  of  the  common  lead- tin  solders  are  around 
400° F ,  this  means  we  cannot  use  ordinary  solder  without  annealing  the 
hard  drawn  copper  tubing.  At  250° F,  the  maximum  allowable  stress  is 
different  for  annealed,  light,  and  hard  drawn  copper  tubing. 

To  be  on  the  safe  side,  we  have  decided  to  specify  that  the 
alloy  to  be  used  to  solder  the  copper  tubing,  in  a  close-packed  cylin¬ 
drical  array,  must  melt  below  200° F.  Table  3  gives  the  composition  of 
some  eutectic  fusible  alloys. 

TABLE  3.  -  Composition  and  -melting  temperatures  of  eutectic  fusible  alloys 


M. 

T. 

Compo 

sition 

Item 

°  F 

°C 

Bi 

Pb 

Sn 

Cd 

Other 

A 

117 

46.8 

44.70 

22.60 

8.30 

5.30 

19.10  In 

B 

,  136 

58  . 

49.00 

18.00 

12.00 

— 

21.00  In 

C 

158 

70 

50.00 

26.70 

13.30 

10.00 

— 

D 

197 

91.5 

51.60 

40.20 

8.20 

. 

* 
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WEIGHT  OF  COPPER  TUBING  IN  INTERCHANGER 


We  have  a  total  of  55  parallel,  3/16-inch  O.D.,  0.042-inch 
wall,  copper  tubes  each  10.5  feet  long.  This  is  a  total  of  577.5 
feet . 

3/16"  =  0.1875  inch 

.1875  -  2  X  .042  =  0.1035 

tt/4  X  (  .  1875)  2  =  0.0276116 
tt/4  X  (  .  1035) 2  =  0.0084134 

The  cross-sectional  area  occupied  by  copper  in  a  single  tube  is 


=  0.0276116  -  0.0084134 
=  0.0191982  square  inches 


Volume 


0.0191982  X  577.5 

144 

=  0.076993  cu.  ft. 

=  133.044  cu.  in. 

=  2180.2  cm3 


Weight  copper  =  2180.2  X  8.5 

=  18,531.7  grams 
=  40.85  lbs. 


APPROXIMATE  OUTSIDE  DIAMETER  OF  PACKED  INTER CHANGER  TUBING 

It  is  planned  to  soft  solder  the  interchanger  tubes  in  a  close- 
packed  cylindrical  array.  For  this  packing,  the  filling  factor  is 


. 


. 
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given  by 

f  =  ^=~  =  0.9069. 

2  2 

The  area  taken  up  per  tube  is  ^ —  =  0.027612 

sq  in. 

There  are  a  total  of  5  +  19  +  31  =  55  tubes,  making  a  total 
area  of 

55  X  .027612  =  1.51866  sq.  in. 

The  interchanger  will  then  take  up  an  area  of 


1.51866 

.9069 


1.6746  sq.  in. 


This  area  is  enclosed  in  a  circle  with  a  diameter  of 


D  =  1.4602  inches. 


ESTIMATION  OF  WEIGHT  OF  LOW  MELTING  SOLDER  IN  INTERCHANGER 

The  cross-sectional  area  of  the  interchanger  is  1.6746  square 
inches.  The  actual  area  taken  up  by  the  copper  tubes  is  1.51866 
square  inches.  The  difference,  0.15594,  is  the  area  available  to  be 
filled  with  solder. 

0.15594  X  10.5 
Volume  of  solder  =  - - - — 


0.01137  cubic  feet 


♦ 


' 
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=  19.647  cubic  inches 

=  321.96  cm^ 


and  with  a  specific  gravity  of  9.58  for  solder,  we  have  a  total 
weight  of 

321.96  X  9.58  =  3084  grams 


or 


3084 

453.6 


6.80  lbs  of  solder. 


ESTIMATION  OF  HEAT  LEAK  AND  WEIGHT  OF  INSULATION 


Linde  insulation  CS-5  is  a  mixture  of  copper  flakes  and  Santocel 
(3}  with  an  apparent  thermal  conductivity  of 


(_3)  Riede,  P.  M.  ,  and  D.  I-J.  Wang,  Characteristics  and  Applications 
of  Some  Super insuiat ion.  Advances  in  Cryogenic  Engineering, 
v.  5,  K„  D.  Timmerhaus,  Editor.  Plenum  Press,  Inc.,  New  York, 
I960,  pp .  209-215. 


-5  -1  -2  -1 

k  =  22  X  10  Btu  hr  ft  (deg  F)  ft  and  a  density  of  11.0  pounds 
per  cubic  foot.  In  order  to  obtain  this  low  value,  it  is  necessary  to 
evacuate  the  insulation.  A  vacuum  of  1  micron  of  mercury  is  suffi- 
c ient . 

We  will  calculate  the  heat  leak  for  a  thickness  of  0.5199  inch. 


The  outside  diameter  of  the  insulation  is  then 


■ 


.q 


' 
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1.4602  +  2  X  .5199  =  2.500  inches 


A. 

i 


tt  X  1.4602  X  10.5 

12 


4.01393 


A 

o 


tt  X  2.5  X  10.5 

12 


6.87223 


A  -  A.  =  2.8583;  A  /A.  =  1.712095 

o  1  o  1 


log  A  /A.  =  0.23353;  In  A  /A.  =  0.53772 

o  1  0  1 


2.8583 

0.53772 


5.3156 


Q 


k  AT  A 
_ av 

L 


We  will  take  the  average  AT  as  being  18° F.  Then 


Q  = 


22  X  IQ"5  X  18  X  5.3156  X  12 

0.5199 


Q  =  0 .486  Btu  hr 


With  20  scfm  of  He  and  and  40  scfm  of  He-^  mixture  going 
through  the  interchanger,  the  heat  exchanged  is  of  the  order  of  1200 
Btu  hr  1  ,  so  the  above  heat  leak  is  less  than  0 .  L°h  of  the  heat  ex¬ 
changed  . 

I 

The  cross-sectional  area  taken  up  by  the  insulation  is  given  by 


tt/4  (2.52  - 


1.46022) 


3.2341  sq.  in. 


18 


The  volume  of  the  insulation  is 


3.2341  X  10.5 
144 


0.2358  cu  ft. 


Then  the  weight  of  the  insulation  is 

0.2358  X  11.0  =  2.59  lbs. 


DETAILS  OF  CALCULATION  OF  TEMPERATURES  IN  INTERCHANGER 
We  have  calculated  the  operating  characteristics  of  this  inter¬ 
changer  under  the  following  operating  conditions. 

Pressure  =  500  psia 

Temperature  of  the  entering  helium  and  nitrogen  streams 
=  20°  C  =  68°  F. 

Temperature  of  the  entering  mixture  stream  =  0^ C  =  32° F. 

Let  stream  1  be  the  helium  stream. 

Let  stream  2  be  the  mixture  stream. 

Let  stream  3  be  the  nitrogen  stream. 

The  fundamental  equation  for  heat  transfer  is  given  by  (4) 


(4)  Walker,  W.  H.,  W.  K.  Lewis,  W.  H.  McAdams  and  E.  R.  Gilliland, 
Principles  of  Chemical  Engineering,  McGraw-Hill  Book  Company, 
Inc.  New  York  and  London,  3rd  ed .  1937,  p.  107. 


' 
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d£ 

dA 


h  (TCu  -  T) 


+  WC 

P 


dT 

dA 


(1) 


where 


d2  , 


dA 


is  the  heat  transferred  to  the  gas  per  unit  area  of  heat  transfer 


-1  -2 

surface  per  unit  time  in  Btu  hr  ft  ,  h  is  the  individual  heat 


transfer  coefficient  for  heat  transferred  from  the  metal  surface  to 

- 1  -2  - 1 
the  gas  in  Btu  hr  ft  (deg  F) 


T^,  is  the  temperature  of  the  copper  tube  surface  in  deg  F. 

T  is  the  temperature  of  the  gas  in  deg  F. 

W  is  the  mass  flow  rate  of  the  gas  in  lbs  hr  ^ . 

is  the  heat  capacity  of  the  gas  in  Btu  lb  ^  (deg  F)  ^ . 

The  +  sign  in  front  of  W  is  to  be  taken  when  the  gas  is  flowing 
in  the  direction  of  dA  positive  and  the  -  sign  is  to  be  taken  when  dA 
positive  is  in  the  opposite  direction  of  gas  flow. 

Now  A  =  A  x  and 

dA  =  A  dx  where  A  is  the  area  available  for  heat  transfer 
per  unit  length  of  interchanger  and  x  is  the  length  of  interchanger. 

In  our  interchanger,  streams  1  and  3  are  flowing  in  the  same 
direction  and  stream  2  is  flowing  in  the  opposite  direction  through 
the  interchanger.  We  then  have  for  the  three  streams 


dQ  dT 

— -  =  A  h.  (T  -  T.)  =  -W.C  ”  ~ 

— 1  i  Cu^  1  1  p^  dx 


dx 


(2) 


, 


' 


20 


^2 

dx 


A2h2(T 


(3) 


A3h3 


(4) 


The  length  of  interChanger  increases  in  the  direction  of  the 
mixture  flow,  stream  2.  We  neglect  the  thermal  lag  through  the  copper 
tube  surface.  Thus,  we  assume  that 


Then 


T  =  T  =  T  =  T 
Cu;[  Cu2  Cu3  Cu 


W.C  ht 
1  P1  ^T1 

V  Cu  V  A  h  dx 


or 


^Cu'V  “  "Vi 


(5) 


(6) 


(7) 


where 


W.C 
1  Pi 


A^i 


(8) 


and 


dx 


Similarly 


a2T2 


(9) 


(10) 


' 


f. 


. 
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(TCu  -  V  "  -  a3T3 


(ii) 


and  eliminating  T  ,  we  have 

U 


(Tx  -  T2)  =  aLT'  +  a2T' 


(12) 


(T3  -  T2)  =  a2T'  +  a3T' 


(13) 


Assuming  there  is  no  heat  leak  into  the  exchanger ,  we  have 


dQ]  dQ  dQ 

3T  +  5r  +  5T  “  0 


(14) 


Therefore,  we  can  write 


W, C  T'  +  W.C  To 

1  Pi  1  3  P3  3 


w.c  i: 

2  p2  2 


(15) 


or 


Vl  +  b3T3 


"  b2T2 


(16) 


where 


b 


(17) 


Equations  (12)  ,  (13)  ,  and  (16)  are  the  differential  equations 


that  fix  the  solution  of  our  heat  transfer  problem.  The  solutions 


' 
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are  (J5) 


(5}  Barieau,  Robert  E.  Design  Equations  for  Multistream  Inter¬ 
changers.  Memorandum  Report  No.  33,  Helium  Research  Center, 
Bureau  of  Mines,  January  1964,  29  pp. 


r-i  +  a„r  _  ctrti. 

[ - ^1  ft  er2X  +  f- - ^  v 

Ll  -  a^^ J  P2  Ll  -  'a1r3J 


1  +  a^r 


T  =  0!  + 

1  2 


sr3x 


(18) 


t2  =  c»2  +  +  Y2er^X 


(19) 


T3  =  “2  + 


-1  +  a2r2- 
.1  -  a3r2 


a  er2x  +  Y  er3x 

P2e  +  Ll  -  a3r3J  V 


(20) 


where  r^  and  r^  are  the  roots  of  the  equation 


(ala3b2  +  a2a3bl  +  ala2b3)r' 


i  b2 (a^  +  a3)  +  b1(a2  -  a3)  +  b3(a2  -  a^) 


+  b2  "  bI  *  b3 


=  0  (21) 


The  constants  a  ,  02,  and  y2  are  fixed  by  the  boundary  condi 
tions  of  the  problem.  The  boundary  conditions  are  fixed  by  the 
temperatures  of  the  streams  when  they  flow  into  the  interchanger, 
our  convention,  stream  2  (mixture)  flows  in  at  x  =  0  and  streams  1 


By 


/ 

' 
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(helium)  and  3  (nitrogen)  flow  in  at  x  =  L,  where  L  is  the  length 
of  interchanger.  Then  we  have  the  three  boundary  conditions 


^2  + 


r  1  “h  aoro-i  t  r' -1-  t 

[rr*3  *  [- 


.1  +  a 2r^ 


air2 


alr3J  2 


y9e 


r3L 


(22) 


cy_  +  0  „  +  v - 

2  p2  '2 


(23) 


L  1  "h  ^2^2~ 

To  =  oi  0  1 

3  2L1  -  a3r2J 


„  r9L 

32e  2  + 


1  +  V31 
a3r3“ 


[- 


Y2e 


r^L 


(24) 


These  equations  are  to  be  used  in  the  following  way:  From  the 
assumed  flow  conditions,  a^,  a^,  a^ ,  b^,  b2  ,  and  b^  are  calculated. 
These  values  are  then  substituted  in  equation  (21)  and  r2  and  r^  are 
determined  as  the  roots  of  equation  (21) .  r^  and  r^  are  then  substi¬ 
tuted  in  equations  (22)  ,  (23)  ,  and  (24)  ,  and  these  equations  are 
solved  for  ,  32?  anc*  Y2  •  Then  these  values  for  ,  32>  and  are 
substituted  in  equations  (18),  (19),  and  (20)  and  T0^ ,  T^ ,  and  T^  are 
calculated.  We  can  then  calculate  the  AT's  at  the  ends  of  the  inter¬ 
changer  . 

Sample  calculation  for  20  scfm  of  helium,  20  scfm  of 

nitrogen,  and  40  scfm  of  helium-nitrogen  mixture 

In  this  section,  I  give  a  sample  calculation  for  20  scfm  of 
helium,  20  scfm  of  nitrogen  and  40  scfm  of  a  50-50  helium-nitrogen 
mixture.  Calculations  were  also  made  for  20  scfm  of  pure  helium 
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only;  20  scfm  of  pure  nitrogen  only;  20  scfm  of  helium,  20/3  scfm 
nitrogen  and  80/3  scfm  of  a  75-25  helium-nitrogen  mixture;  and  for 
20/3  scfm  helium,  20  scfm  nitrogen,  and  80/3  scfm  of  a  25-75  helium- 
nitrogen  mixture.  Calculations,  other  than  for  the  50-50  mix,  are 
summarized  in  the  tables. 

In  my  calculations,  I  have  taken  C  =  1.23960  Btu  lb  ^(deg  F) 

Pi 

for  helium,  and  C  =  0.26274  for  nitrogen.  For  the  mixture,  I  have 

P3 

assumed,  at  500  psia, 


W_C 
2  P, 


W.C  +  W_C 
1  Pl  3  p3 


=  1  -  0 . 04y ( 1-y) 


(25) 


where  y  is  the  mole  fraction  of  helium  in  the  mixture. 
Then  for  the  50-50  mix 


W_C 
2  p. 


W.C  +  W0C 
1  PL  3  P3 


0.99 


(26) 


I  take  a  standard  cubic  foot  as  the  amount  of  gas,  at  1  atmos¬ 
phere  and  70° F,  in  a  cubic  foot 


1  cu.  ft.  =  28,316.8  cm' 


(27) 


The  molal  volume  of  a  standard  cubic  foot  of  gas  at  70  F  and 


1  atmosphere  is 


294  26  3 

=  22,414.6  X  =  24,146.9  cm 


273.15 


(28) 


. 
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1  SCF 


28.316.8 

24.146.9 


1 . 17269  g  moles . 


M.  W.  He  =  4.0028 
M.  W.  N2  =  28-016 


1  SCFM 


1.17269  X  60  X  M.W. 

453.592 


lbs  hr 


1  SCFM  =  0.155120  X  M.  W.  lbs  hr 

20  SCFM  (He)  =  12.4184  lbs  hr'1  =  W 

20  SCFM  (N  )  =  86.9170  lbs  hr'1  =  W 

40  SCFM  (50-50,  He  -  N  )  =  99.3354  lbs  hr'1  =  W2 


(29) 

(30) 

(31) 


Then 


b.  =  W.C  =  15.3938  Btu  hr'1  (deg  F)'1 

1  1  P1 


(32) 


b  =  W.C  =  22.8366  Btu  hr'1  (deg  F)'1 

3  3  p3 


(33) 


W2Cp  =  b2  =  0.99  (bL  +  b3)  =  37.8481  Btu  hr'1  (deg  F) " 1  (34) 


gg'3331  =  0.381013  Btu  lb'1  (deg  F)'1 


(35) 


Table  4  gives  data  on  the  various  flow  rates  for  the  different 


calculations . 


■  )2  j 


’ 


. 
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TABLE  4.  -  Flow  rates  and  heat  capacities  for  various  assumed 

flow  conditions 


Stream 

Ga  s 

SCFM 

W,  lbs.  hr"1 

Cp,Btu  lb" i(deg  F)  ~ 1 

1 

He 

20 

12.4184 

1.23960 

2 

He ( 100) 

20 

12.4184 

1.23960 

3 

— 

0 

i 

He 

20 

12.4184 

1.23960 

2 

He-ISL  (75-25) 

80/3 

41.3907 

0.551658 

.  3 

N2 

20/3 

28.9723 

0.26274 

1 

He 

20 

12.4184 

1.23960 

2 

He-N_ (50-50) 

40 

99.3354 

0.381013 

3 

N2  2 

20 

86.9170 

0.26274 

1 

He 

20/3 

4.1395 

1.23960 

2 

He-N  (25-73) 

80/3 

91.0565 

0.304845 

3 

N2 

20 

86.9170 

0.26274 

i 

_  __ 

0 

2 

N  ( 100) 

20 

86.9170 

0.26274 

3 

N2 

20 

86.9170 

0.26274 

■ 


, 


, 
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Heat  Transfer  Coefficients  for  Gases 

in  Circular  Tubes 

The  heat  transfer  coefficient,  for  a  gas  under  turbulent  flow, 
inside  a  circular  tube  is  given  by  (6) 


(6)  Reference  4 p .  112. 


(36) 


where 


-1-2  - 1 

h  is  the  heat  transfer  coefficient  in  Btu  hr  ft  (deg  F) 

-1  -2  -1 

k  is  the  thermal  conductivity  in  Btu  hr  ft  (deg  F)  ft. 


D  is  the  inside  diameter  of  the  tube  in  feet 

G  is  the  mass  velocity  ip  lbs  hr  ^  (sq.  ft.  of  cross-section)  ^ 

-1  -1 

is  the  viscosity  in  lbs  hr  ft  =  2.42  X  centipoises 
Cp  is  heat  capacity  in  Btp  lb  ^  (deg  F)  ^ 


Walker,  Lewis,  McAdams,  apd  Gilliland  (7_)  say  that  for  the 

(_7)  Reference  4,  p.  113. 

common  gases  equation  (36)  can  be  simplified  to  give  the  dimensional 
equation 


G 


0.8 


h  =  0.0144  Cp 


0.2 


(37) 


D 


. 
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However,  Giauque  ( 8 )  obtained  the  expression 


(J3)  Giauque,  W.  F.  Liquid  Oxygen  Trailer  Unit,  Design,  Construction, 
and  Operation.  Final  Report  to  National  Defense  Research  Com¬ 
mittee.  Report  OSRD  No.  4141,  July  1944,  p.  V-ll>  539  pp. 


0.8 

h  =  0.0120  Cp  (38) 

DU*Z 


for  high  pressure  air,  inside  straight  circular  tubes.  We  have 
accepted  equation  (38)  for  calculating  ail  heat  transfer  coefficients. 


G 


W4 

2 

nrrD 


(39) 


where  n  is  the  number  of  parallel  tubes  and 


.0.8  TT0 . 8 

4 _ W _ 

0.8  0.8  1 . 6 

n  tt  D 


(40) 


so  that 


,0.8  0.8 

h  =  0.0120  Cp  -5^  (  -  )  - 

tt  D 


4 

TT 


0.8 

O 


1.213188 


Thus 


W  0-8 

h  =  0.0145583  (  -  ) 

n 


C£ 


(41) 


(42) 


(43) 


■ 
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For  3/16"  O.D.  tubing  with  0.042"  wall 


I.D.  =  .1875  - 

.084  =  0.1035  inch  (44) 

0.1035 

12 

=  0.008625  ft.  (45) 

i  -  115- 

,942029  ft"1  (46) 

ii 

00 

Q 

5195.575  (47) 

Then 

h  =  75.6387 

W  ^  ^ 

(  l  )  Cp  (48) 

For  helium,  n^  =  5 

For  nitrogen,  n^  =  19 

For  mixture,  n  =31 

2' 


nL°  =  3.623898 

Cp^  =/  1.2396 

0.8 

n2 

=  15.598734 

n3°'8  =  10.543939 

Cp3  =  0.26274 

Then  we  have  for  the  helium  at 

P  =  500  psia, 

h  =  25.8732  W^'8 


(49) 


' 
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and  for  the  mixture  at  P  =  500  psia, 


ii 

CNJ 

rC 

4.84903  Cp2  W20,8 

(50) 

and  for 

the  nitrogen  at  P  = 

:  500  psia, 

h3 

1.884809  W3°-8 

(51) 

With  W 

=  12.4184,  W2  =  99. 

3354,  W3  =  86.9170,  and  Cp2 

=  0.381013, 

we  have 

w  °-8 
wi 

=  7.5033035 

(52) 

„  0.8 
W2 

=  39.598910 

(53) 

W  °'8 

=  35.586378 

(54) 

Finally,  we  have 


hL  =  194.1345  Btu  hr'1 

ft  2  (deg  F)  1 

(55) 

h2  =  73.1607  Btu  hr'1 

ft'2  (deg  F)'1 

(56) 

h  =  67.0735  Btu  hr”1 

ft'2  (deg  F)'1. 

(57) 

Now 

WCp  =  WCp 
a  Ah  nrrDh 

ttD  =  tt  X  0.008625  =  0.0270962366 


. 
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36.905494  WCp 


Stream 

n 

h 

WCp  =  b 

nh 

a 

1 

5 

194.1345 

15.3938 

970.673 

0.585280 

2 

31 

73.1607 

37.8481 

2267.982 

0.615879 

3 

19 

67.0735 

22.8366 

1274.397 

0.661329 

The 

var ious 

values  of  h, 

b ,  and  a  are 

summar ized 

in  table  5. 

then  have 


a  =  0.585280 

b.  = 
i 

15.3938 

a2  =  0.615879 

b2  ' 

37.8481 

a3  =  0.661329 

b3  " 

22.8366 

a1a3b2 

=  14.6495854 

a2a3bl 

=  6.26987385 

aia2b3 

=  8.23171877 

ala3b2  +  a2a3bl 

+  axa2b3  =  29 

.15117802 

(58) 

b 

b 

b 


2 

1 

3 


(aL  + 

a3^ 

=  47.18178209 

(a2  " 

a3) 

=  -0.69964821 

(a2  " 

al} 

=  0.69877712 

b2(ai  +  a3)  +  bi  (a: 


-  a3)  +  b3  (a2  -  aL)  =  47.18091100 


(59) 


’ 
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TABLE  5. 

-  Various 

values  of  h,  b,  and  a  for  the 

various  flow 

Stream 

SCFM  h,Btu  hr  ‘'"ft  ^  (deg  F) 

-1  b 

a 

1 

20 

194.1345 

15.3938 

0.585280 

2 

20 

45.1013 

15.3938 

0.406337 

3 

0 

1 

20 

194.1345 

15.3938 

0.585280 

2 

80/3 

52.5832 

22.8335 

0.516957 

3 

20/3 

27.8518 

7.6122 

0.530878 

i 

20 

194.1345 

15.3938 

0 , 585280 

2 

40 

73.1607 

37.8481 

0.615879 

3 

20 

67.0735 

22.8366 

0.661329 

1 

20/3 

80.6136 

5.1313 

0.469829 

2 

80/3 

54.5988 

27.7581 

0.605252 

3 

20 

67.0735 

22.8366 

0.661329 

1 

0 

2 

20 

45.3383 

22.8366 

0.599647 

3 

20 

67.0735 

22.8366 

0.661329 

• 

1 

' 
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b2  -  bi  -  b3  =  -0.3823 


or 


Substituting  in  equation  (21)  ,  we  have 

29.15117802  r2  -  47.180911r  -  0.3823  =  0 


r2  -  1.6184907164  r  -  0.013114392829  =  0 


The  roots  of  this  equation  are 


r2  =  1.626553404 

r3  =  -0.0080626880 


(60) 

(61) 

(62) 


The  equations  for  r  and  the  roots  for  the  various  calcula¬ 
tions  carried  out  are  given  in  table  6. 


TABLE  6. 

-  Equation 

for 

r 

and  roots  for  the  various  calculations 

carried  out 

Stream 

SCFM 

1 

20 

2 

r 

— 

1.8169475284  r 

-  0.012662853174  =  0 

2 

80/3 

r? 

= 

1.823890300 

3 

20/3 

z 

r3 

= 

-0.0069427715 

1 

20 

2 

r 

_ 

1.6184907164  r 

-  0.013114392829  =  0 

2 

40 

r? 

= 

1.626553404 

3 

20 

z 

r3 

= 

-0.0080626880 

1 

20/3 

2 

r 

_ 

1.991757261  r  - 

0.01221713692  =  0 

2 

80/3 

ro 

= 

1.997872335 

3 

20 

z 

r3 

= 

-0.006115074 

' 

' 


6  j  d&j  ivW?  sii  3|#o  b-jty  <*.nol  i 
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Then 


a^2  =  0.95198917629;  1  -  a  r  2  =  0.04801082371 


a2r  2  =  1-0017600839;  1  +  a2r  =  2.0017600839 


a3r2  =  1.07568693611;  1  -  a^  =  -0.07568693611 


1  +  a2r2 
1  -  a!r2 


41.69393335 


1  +  a2r2 
1  "  a3r2 


-  26.44789427 


a][r3  =  -0.004718930033;  1  -  a  r  =  1.00471893003 


a2r3  =  -0.004965640223;  1  +  a  r3  =  0.995034359777 


,rn  =  -0.005332089392;  1  -  anr0  =  1.005332089392 


a3"  3  - *  ~r3 


1  +  a  r  1  +  a  r 

Z  J  =  0.99036091593;  - - —  =  0.98975688796 


1  -  aLr3 


1  -  a3r3 


Substituting  in  equations  (18)  ,  (19)  ,  and  (20)  ,  we  have 


T  =  a2  +  41.69393335  3^2*  +  0 . 9903609 159372eL (63) 


r  ox 


r0x  r3X 


T2  =  a2  +  P2e‘2"  +  Y2e 


(64) 


J 
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T3  =  °^2  "  26-44789427  ^2^^  +  °-98975688796v2er3X  (65) 


With  the  length  of  the  interchanger  being  10.50  feet,  we  have 


r2L 


r£L  =  17.078810742;  2  303  =  7-4172333 


(66) 


r3L 


r3L  =  -0.084658224;  =  -0.0367666 


(67) 


[10] 


v_ 

2.303 


r„L 


=  e 


2  =  2.6135650  X  10? 


(68) 


[10] 


"r3L 

2.303 


e"r3L  =  1.0883450 


(69) 


r  „L 

e  3  =  0.91882629129 


(70) 


We  then  have  the  three  boundary  condition  equations, 


68  =  c*2  +  41.69393335  X  2.613565  X  10  $2 


(71) 


+  0.99036091593  X  0.91882629129  y. 


32  =  a2  +  B2  +  Y2 


(72) 


or 


68  =  a2  -  26.44789427  X  2.613565  X  10  32 


+  0.98975688796  X  0.91882629129  Y, 


(73) 


=  c*_  +  1.0896980492  x  10y  39  +  0.9099696474  y  (74) 


68 


■ 
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32  =  c*2  +  g2  +  y2  (75) 

68  =  c*2  -  0.6912329079  X  109  P2  +  0.90941465069  y2  (76) 

Subtracting  equation  (76)  from  (74)  ,  we  have 
0  =  1.780930957  X  i0932  +  0 . 0005549967y2  (77) 

and  subtracting  equation  (75)  from  (74)  ,  we  have 

36  =  1.0896980482  X  109  $2  -  0.0900303526  y2  (78) 

From  equation  (77) 

-13  , 

0  =  02  +  3.116329118  X  10  Y2  (79) 

and  from  equation  (78) 

33.036674761  X  10“9  =  -  826.19541041  X  10' I3Y2  (80) 

Subtracting  equations  (79)  and  (80) 

33.036674761  X  10“9  =  -829.3117395  X  10“13y2  (81) 

Y2  =  -398.36256  (82) 

P2  =  1.241429  X  i0”10  (83) 


' 


' 
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and 

c*2  =  430.36256 

Then  we  have 

T  =  430.3626  +  5.1760  X  lO'V^  -  394.5227  e^  (84) 

T2  =  430.3626  +  1.2414  X  10_10er2X  -  398.3626  e^  (85) 

T3  =  430.3626  -  3.2833  X  10"9er2X  -  394.2821  er3X  (86) 

Equations  (84)  ,  (85) ,  and  (86)  give  the  temperatures  of  the 
streams  through  the  interchanger.  When  x  =  L  =  10 . 50  feet,  we  have 

=  430.3626  +  0.1353  -  362.4978  =  68.0000°F  (87) 

T2  =  430.3626  +  0.0032  -  366.0260  =  64.3398°F  (88) 

=  430.3626  -  0.0858  -  362.2768  =  68.0000°F  (89) 

At  the  cold  end  of  the  interchanger ,  where  x  =  0 ,  we  have 
’f  =  430.3626  -  394.5227  =  35. 8399° F  (90) 

t2  =  430.3626  -  398.3626  =  32.0000°F  (91) 

=  430.3626  -  394.2821  =  36. 0805° F  (92) 


' 
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Therefore , 


=  3 . 660°  F 


(93) 


or 


2 . 034°  C 


temperature  of  the  entering  pure  helium  and  nitrogen.  The  helium 
stream  will  be  cooled  to  35.840° F  or  to  3.840°F  of  the  tempera¬ 
ture  of  the  entering  cold  mixture  stream. 

While  the  helium  and  nitrogen  streams  start  out  at  the  same 
temperature  at  the  warm  end  of  the  interchanger,  when  they  come 
out  of  the  interchanger  at  the  cold  end,  the  helium  stream  is 
0.241° F  colder  than  the  nitrogen  stream.  The  problem  of  allowing 
or  correcting  for  this  temperature  difference  will  be  considered 
later.  The  temperatures  of  the  various  streams  as  a  function  of 
interchanger  length  are  given  in  table  7.  The  temperature  of  the 
copper  tube  wall  surface  as  a  function  of  interchanger  length  is 
also  given. 

Table  8  gives  the  temperatures  calculated  at  the  hot  and 
cold  end  of  the  interchanger  from  the  equations  given  in  table  7. 
The  last  column  in  table  8  gives  the  temperature  difference  be¬ 
tween  the  streams  at  the  warm  end  of  the  interchanger.  This 
temperature  difference  would  be  zero  in  an  infinite  interchanger. 


Table  9  gives  the  rate  of  heat  transfer  in  the  inter changer . 
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TABLE 

Stream 

1 

2 


1 

2 

3 


1 

2 
3 


1 

2 

3 


2 

3 


.  -  Temperature  of  streams  and  tube  wall  surface  as  a  function 

of  interchanger  length 


SCFM 

20  T  =  35.1065  +  3.13272x 

20  T2  =  32.0000  +  3.13272x 

Tn  =  33.2729  +  3.13272x 
Cu 

20  T  =  497.8375  -  2.4925  X  10“10er2X-  462.2870er3X 

80/3  T2  =  497.8375  +  0.0866  X  10~i0er2x-  465.8375er3X 

20/3  T3  =  497.8375  +  5.3003  X  10_10er2X-  462.46i0er3X 

T  =  497.8375  +  0.1682  X  10"10er2X-  464.1655er3X 
Cu 

20  T  =  430.3626  +  51.760  X  10"10er2X-  394.5227er3X 

40  T2  =  430.3626  +  1.2414  X  10'10er2X-  398.3626er3X 

20  T3  =  430.3626  -  32.833  X  10~i0er2X-  394.2821er3X 

T_  =  430.3626  +  2.4850  X  10"10er2X-  396.3844er3X 
Cu 

20/3  T  =  551.85215  +  36.777  X  iO_11er2X-  516 .44433er3X 

80/3  T2  =  551.85215  +  1.021  X  10_11er2X-  519 . 85215er3X 

20  T3  =  551.85215  -  7.023  X  lO^V2^  515 . 84199er  3X 

T^  =  551.85215  +  2.256  X  10“Uer2X-  517 . 92809er3X 
Cu 

20  T2  =  32.0000  +  3.06097x 

20  T3  =  35.8598  +  3.06097x 

T  =  33.8355  +  3.06097x 
Cu 


' 


- 

/ 


. 
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TABLE  8 .  -  Temperatures  at  hot  and  cold  end  of  interchanger 


Stream 

SCFM 

T° 

tl 

TL  -  T° 

TL  -  T 
1  : 

1 

20 

35.1065 

68.0000 

32.8935 

3.1065 

2 

Cu 

20 

32.0000 

33.2729 

64.8935 

66.1665 

32.8935 

1 

20 

35.5505 

68.0000 

32.4495 

3.2473 

2 

80/3 

32.0000 

64.7527 

32.7527 

3 

Cu 

20/3 

35.3765 

33.6720 

68.0000 

66.3089 

32.6235 

1 

20 

35.8399 

68.0000 

32.1601 

3.6602 

2 

40 

32.0000 

64.3398 

32.3398 

3 

Cu 

20 

36.0805 

33.9782 

68.0000 

66.1607 

31.9195 

1 

20/3 

35.4078 

68.0000 

32.5922 

3.6570 

2 

80/3 

32.0000 

64.3430 

32.3430 

3 

Cu 

20 

36.0102 

33.9241 

68.0000 

66.1633 

31.9898 

2 

20 

32.0000 

64.1402 

32.1402 

3.8598 

3 

Cu 

20 

35.8598 

33.8355 

68.0000 

65.9757 

32.1402 

. 


' 

■ 

TABLE  9. 


Calculation  of  heat  transferred 


Stream 

SCFM 

TL  -  T° 

b 

Q,  Btu  hr 

1 

20 

32.8935 

15.3938 

-506.356 

2 

20 

32.8935 

15.3938 

+506.356 

1 

20 

32.4495 

15.3938 

-499.521 

2 

80/3 

32.7527 

22.8335 

+747.859 

3 

20/3 

32.6235 

7.6122 

-248.337 

Z  =  +0.001 

1 

20 

32.1601 

15.3938 

-495.066 

2 

40 

32.3398 

37.8481 

+1224.000 

3 

20 

31.9195 

22.8366 

-728.933 

Z  =  +0.001 

1 

20/3 

32.5922 

5.1313 

-167.240 

2 

80/3 

32.3430 

27.7581 

+897.780 

3 

20 

31.9898 

22.8366 

-730.538 

Z  =  +0.002 

2 

20 

32.1402 

22.8366 

+733.973 

3 

20 

32.1402 

22.8366 

-733.973 

, 

.. 

. 

. 
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CALCULATION  OF  TEMPERATURE  DROP  THROUGH  COPPER  TUBE  WALL 

The  equations  used  in  calculating  the  performance  of  the 
interchanger  have  been  derived  on  the  assumption  that  the  tempera¬ 
ture  drop  through  the  copper  wail  is  small  and  can  be  neglected. 

It  is  important  that  this  temperature  difference  be  calculated  to 
see  if  it  is  permissible  to  neglect  it. 

From  table  9,  we  see  that  for  the  50-50  mix  experiment, 


Q  =  -495.066  Btu/hr 


and  since  for  the  helium  stream  there  are  five  tubes,  we  have 


99.013  Btu/hr  tube. 


The  area  inside  a  tube  is 


ttD.L 


A. 

l 


l 


where  D  is  the  inside  dianieter  in  feet  and  L  is  the  length  of 
tube  in  feet 


A 


it  X  0.008625  X  10.50 


l 


A.  =  0.28451  ft 


2 


l 


y  g  |<j 

, 

* 
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The  area  outside  a  tube  is 


A  =  ttD  L 
o  o 


it  X  .1875 

12 


X 


10.50 


A  =  0.51542  ft2 

o 

A  -  A. 

A  =  -2 - i_ 

av  In  A  /A. 

o  1 


A  -  A.  =  0.23091 

o  l 

A  /A.  =  1.811606 

o  1 

iog.„  A  /A.  =  0.25806 

10  o  l 


In  A  /A.  =  0.59421 

o  1 

A  =  0.38860 


av 


Q 


k  A  AT 
av 


Cu 


L 


.042 

12 


0.0035  ft. 


_£L_ 

kA 

av 


ATcu  =  Q  x 


0.0035 


222  X  0.38860 


AT  =  Q  X  0.0000406 
Cu 


X1  - 1 

Then  for  the  helium,  with  —  =  -99.013  Btu  hr  tube 

ni 


we  have 


AT  =  0 . 0040  F 

Cui 


Table  9  indicates  that  for  the  nitrogen  stream 


Q3  =  -728.933  Btu  hr 


-1 


and  with  19  parallel  tubes  for  the  nitrogen  stream,  we  have 


-  38.365  Btu  hr  ^  tube  ^ 


n. 


and  thus,  for  the  nitrogen  stream 


AT_,  =  38.365  X  .0000406 

Cu3 


AT  =  0 . 0016  F 

Cu3 


Table  9  indicates  that  for  the  50-50  mixture  stream 


=  1224.000  Btu  hr 


-1 
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and  with  31  parallel  tubes  for  the  mixture  stream,  we  have 


39.484  Btu  hr  ^  tube 


Thus,  for  the  mixture  stream 


AT 


39.484  x  0.0000406 


Cu 


2 


AT 


Cu 


0 . 0016°  F 


2 


With  the  AT  in  the  interchanger  between  the  streams  being 

cooled  and  heated  being  3.7°F,  we  see  that  the  AT  through  the 

Cu 


copper  tube  wall  surface  is  completely  negligible. 

DETAILS  OF  CALCULATION  OF  PRESSURE  DROPS  THROUGH  THE  INTERCHANGER 
If  we  neglect  kinetic  energy  and  gravitational  potential  energy 
effects,  the  pressure  drop  due  to  friction  for  circular  tubes  can  be 
expressed  by  the  so-called  Fanning  equation  (_9) 

(9)  Reference  1,  p.  77. 


dF 

dN 


dp  JL 
dN  p 


(94) 


. 


' 


' 

' 
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where  p  is  the  pressure  in  pounds  per  square  foot 
p  is  the  density  in  pounds  per  cubic  foot 
N  is  the  length  in  feet 
f  is  the  dimensionless  friction  factor 
V  is  the  average  linear  velocity  in  feet  per  second 
g  is  the  acceleration  due  to  gravity 
g  =  32.17405  ft  sec  ^ 

D  is  the  inside  diameter  of  the  circular  pipe  in  feet. 


For  gas  flow,  we  replace  p  by  where  v  is  the  spe 

RTZ 

in  cubic  feet  per  pound,  v  is  then  replaced  by  v  =  , 

the  universal  gas  constant  in  f t- lb/lb-mole  deg  F. 


cific  volume 
where  R  is 


R  =  1545.393  ft- lb/ lb-mole  deg  F 
T  is  the  absolute  temperature  in  degrees  Rankine 
Z  is  the  dimensionless  compressibility  factor 
M  is  the  molecular  weight  of  the  gas 
p  is  the  pressure  in  pounds  per  square  foot 

V  =  (95) 

where  w  is  the  total  mass  rate  of  flow  in  lbs  per  second 
n  is  the  number  of  parallel  tubes 

S  is  the  cross-sectional  area  at  right  angles  to  the  tube  in 


square  feet  per  tube 
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S  = 


rrD 

4 


(96) 


V 


w 


2  2  Q  2 

n  p  S 


(97) 


2n4 

TT  D 

16 


(98) 


V 


16 _ w_ 

2  4  2  2 

tt  D  p  n 


(99) 


Substituting  in  equation  (94) ,  we  have 


_!  d_p 

f  dN 


32  w  p 


2  4  22 
tt  D  p  n  gD 


32  w2y 

2  n5  2 
tt  gD  n 


(100) 


and  substituting  for  v,  we  have 


1  d£  =  32RTZw 

f  dN  2Mn52 

TT  gMpD  n 


(101) 


32R 

2 

TT  g 


32  X  1545.393 


(3.14159)  X  32.17405 


=  155.734263 


(102) 


Then  for  gas  flow  in  general, 


_1  d_p 
f  dN 


155.734263  TZw2 

MpD5n2 


(103) 
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In  our  interchanger , 


3/16  -  2  X  .042 

12 

D  =  0.008625  ft.  =  0.8625  X  10_2ft. 


0.477304  X  10 


-10 


Then 


1_  dj) 
f  dN 


326.279  X  1010  TZw2 


Mpn' 


(104) 


As  our  lowest  operating  pressure  will  be  500  psia,  we  will  take 
p  =  500  X  144  =  72000  psfta, 

and  we  will  take 

T  =  283.15  X  1.8  =  509 . 67° R 


so  that 

326.279  X  101Q  X  509.67 

72000 


2.30965  X  1010 


Thus  , 


1  dp 
f  dN 


10  2 

2.30965  X  10  Zw 

Mn2 


(105) 


For  stream  3,  the  nitrogen  stream,  Z^  at  500  psia  =  0.99149 


=  28.016  and  n^  = 


19 


. 


' 
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361 


2.30965  X  1010  X  0.99149  w^ 
28.016  X  361 


% 

dN 


2.26423  X  106  w^ 


(106) 


(107) 


For  stream  l,  the  helium  stream 


Z  =  1.01639 


M  =  4.0028 


=  25 


so  that 


1  dp  2.30965  X  1010  X  1.01639  w^ 

f~  dN  “  4.0028  X  25 


^1 

dN 


2.34586  X  108 


(108) 


(109) 


For  stream  2,  the  mixture  stream,  with  50-50  mix,  we  will  take 

z2  =  1.0000 


28.016  +  4.0028 

2 


16.0094 


31 


2 

ru 


961 


50 


Thus  , 


L  dp2  2.30965  X  1010  w2 
T  dT  =  16.0094  X  961 


(110) 


1.50123  X  106  w2 


(111) 


In  order  to  apply  equations  (107)  ,  (109)  ,  and  (111)  ,  it  is 
necessary  to  read  the  friction  factor,  f,  from  a  graph  given  on  page 
78,  of  reference  1.  In  order  to  do  this,  it  is  necessary  to  know  the 
Reynold  1 s  number ,  Re , 


Re 


4  w 

rrD|in 


(112) 


where  w  is  the  mass  flow  rate  in  lbs  per  sec.  and  D  is  the  I.D.  of  the 
tube  in  ft. 


(j,  is  the  viscosity  in  Ib/sec  ft. 

p,  =  0.0672  times  the  viscosity  in  poises 

n  is  the  number  of  parallel  tubes. 

D  =  0.8625  X  10~2  ft. 


4_ 

ttD 


147.622 


Therefore , 


147.622  w 

(in 


Re 


(113) 
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I  have  taken  the  viscosity  of  helium  and  nitrogen  from 
Richardson,  Gordon,  Cooper ,  and  Walker  (10)  as  being 

( 10)  Richardson,  H.  P.,  J.  L.  Gordon,  J.  L.  Cooper  and  J.  D.  Walker. 
Thermophysical  Properties  of  Selected  Gases  Below  300°K.  He¬ 
lium  Research  Center  Internal  Report  No.  34,  July  1963. 


T)^e  =  192  micropoises  at  500  psia  and  10  C 


T)n  =  177  micropoises  at  500  psia  and  10  C 

uu  =  192  X  10’6  X  0.0672  -  1.29  X  10"5 
pHe 

*  177  X  10’6  X  0.0672  -  1.19  X  10'5 


^2 *  *  ^ (He  -  N2) 


1.24  X  10 


-5 


For  the  50-50  mixture,  I  take  the  mean  of  the  pure  component  values 


Then  with  n^  =  5 , 


147.622  w. 


Re 


1  1.29  X  10"5  X  5 


Re ^  =  2.2887  X  10  wx 


(114) 


f 


■ 
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With  n0  =  31 , 

147.622  w£ 

1.24  X  10"5  X  31 

3.8403  X  105  w2  (115) 

147.622  w3 
1.19  X  10”5  X  19 

6.5291  X  105w3  (116) 

I  summarize  pressure  drop  calculations  for  various  flow  rates  in 
tables  10  and  11  of  this  report. 

ESTIMATION  OF  ADDITIONAL  INTERCHANGER  TUBING  NEEDED  TO  HAVE  THE  COPPER 

TUBE  WALL  INDICATE  THE  FINAL  EQUILIBRIUM  TEMPERATURE  OF  THE  HELIUM- 

NITROGEN  MIXTURE  STREAM  TO  WITHIN  0.016F~ 

At  the  warm  end  of  the  interchanger,  the  copper  tube  wall  surface 
will  be  at  a  temperature  between  that  of  the  entering  gases  and  the 
exiting  mixture  gas  stream.  We  ask  the  question,  "How  much  additional 
length  of  tubing  is  required  to  have  the  temperature  of  the  tube  wall 
within  0.01° F  of  the  mixed  gas  temperature?" 

In  order  to  answer  this  question,  it  is  necessary  to  allow  for 
heat  conduction  along  the  tube  wall.  If  x  represents  the  variable 
length  of  the  interchanger  ,  we  assume  that  the  temperature  of  the  tube 


Re, 


Re, 


and  with  n3  =  19, 


Re. 


Re. 


' 
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TABLE 

10.  - 

Pressure  drop  calculations  for  various  flow 

rates 

Stream 

SCFM 

w,  lbs  sec  ^ 

Re 

f 

_ 

dN 

1 

20 

3.44956  X  10"3 

7,895 

0.0082 

22.89 

2 

20 

3.44956  X  10"3 

1,273 

0.013 

0.944 

1 

20 

3.44956  X  iO"3 

7,895 

0.0082 

22.89 

2 

80/3 

11.49742  X  10"3 

4,328 

0.0098 

3.112 

3 

20/3 

8.04786  X  10'3 

5,255 

0.0092 

1.349 

1 

20 

3.44956  X  10"3 

7,895 

0.0082 

22.89 

2 

40 

27.59316  X  10“3 

10,597 

0.0077 

8.801 

3 

20 

24.14361  X  10"3 

15,764 

0.0069 

9.107 

1 

20/3 

1.14986  X  10‘3 

2,632 

0.0125 

3.877 

2 

80/3 

25.29347  X  10'3  ♦ 

9,913 

0.0078 

5.448 

3 

20 

24.14361  X  10-3 

15,764 

0.0069 

9.107 

2 

20 

24.14361  X  10"3 

9,661 

0.0078 

3.867 

3 

20 

24.14361  X  10"3 

15,764 

0.0069 

9.107 

TABLE 

11.  - 

Pressure  drop  calculations  for  various  flow 

rates , 

N  =  10.5 

ft 

Stream 

SCFM 

Ap,  psf. 

Ap ,  psi 

1 

20 

240.3 

1.669 

2 

20 

9.91 

0.069* 

* Laminar  Flow 

1 

20 

240.3 

1.669 

2 

80/3 

32.68 

0.227 

3 

20/3 

14.16 

0.098 

1 

20 

240.3 

1.669 

2 

40 

92.41 

0.642 

3 

20 

95.62 

0.664 

1 

20/3 

40.71 

0.283 

2 

80/3 

57.20 

0.397 

3 

20 

95.62 

0.664 

2 

20 

40.61 

0.282 

3 

20 

95.62 

0.664 

1 

* 

• 
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wall  is  only  a  function  of  x.  We  thus  neglect  radial  temperature 
gradients.  In  a  differential  length  of  interchanger,  we  have 


dQ2  dT 

-  =  n  ttD  h9  (T  -  TJ  =  W  C  — ^ 

2  2  2  Cu  2  2  p2  dx 


dx 


(117) 


TCu  '  T2  "  a2T2 


(118) 


where 


W9C 

2  p. 


a2  n2nD2h2 


(119) 


The  heat  flowing  into  the  differential  length,  dx ,  is  given  by 


dT 


Q .  =  -  kA 


Cu 


in 


dx 


(120) 


where  k  is  the  thermal  conductivity  of  copper  and  A  is  the  cross- 
sectional  area  available  for  heat  conduction.  The  heat  flowing  out 
of  this  differential  length  is  given  by 


dQ, 


‘out 


Q.  + 

in 


m 


dx 


dx 


(121) 


dT 


=  -  kA  ( 


Cu 


d2T 


dx 


+ 


Cu 


dx) 


(122) 


dx 


The  heat  transferred  to  the  gas  is  given  by 


' 
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or 


Q.  -  Q 

in  out 


=  dfi 

dx 


d2T 

dx  =  kA  - ^  dx 

dx 


(123) 


or 


d£ 

dx 


2 

d  T 


=  kA 


Cu 


dx" 


(124) 


dx 


dQ; 

dx 


(125) 


and 


d2T 


kA 


Cu 


dx' 


dQ; 

dx 


dT, 


=  W,C 


2  p2  dx 


(126) 


kAT"cu  ■  b2T2 


(127) 


We  try,  as  solutions  to  equations  (118)  and  (127)  , 


^2  °^2e 


rx 


(128) 


Tr  =  «  e 
Cu  Cu 


rx 


(129) 


m  /  rx 

T2  =  Q^re 


(130) 


' 
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TCu  “  T2  ^Cu  "  a2*e 


rx 


(131) 


T"  =  a  r2erx 
Cu  Cu 


(132) 


Substituting  in  equation  (118)  ,  we  have 


IT  X  IT  X 

(“Cu  ‘  a2)e  =  a2Ve 


(133) 


and  for  this  equation  to  be  true,  it  is  necessary  and  sufficient  that 


“cu  a2  (‘1  +  a2r‘) 


(134) 


Substituting  in  equation  (127)  ,  we  have 


.  .  2  rx  rx 

kA  are  =  b_or0re 
Cu  2  2 


(135) 


and  for  this  equation  to  be  true,  it  is  necessary  and  sufficient  that 


kA  &  r  -  b0a0 
Cu  2  2. 


=  0 


(136) 


Substituting  from  equation  (134) ,  we  have 


kAr  (1  +  a^r)  -  b 


=  0 


(137) 


This  is  a  cubic  equation  with 


0. 


(138) 


- 

' 


' 
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The  other  two  roots  are  given  by 


kAa^r  +  kAr  -  =  0. 


(139) 


There  are  31  parallel  tubes  carrying  the  mixed  stream;  therefore, 
the  cross-sectional  area  for  heat  conduction  is 


A  =  -1—X—Q =  0.00413296148  ft2 
144 


For  copper 


k  =  222  Btu  hr"1  ft“2(deg  F)_1ft 


kA  =  0.917517 


For  the  50-50  mixture  stream, 

a2  =  0.615879 

b2  =  37.8481 

kAa2  =  0.5650794524 

Substituting  in  equation  (139),  we  have 

0.5650794524r2  +  0.917517r  -  37.8481  =  0  (140) 
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or 


r2  +  1 . 623695564r  -  66.9783689  =  0 


(141) 


The  roots  of  this  equation  are 


r2  =  -9.036047785 


(142) 


r3  =  +7.412352225 


(143) 


Only  negative  or  zero  roots  will  fit  our  boundary  conditions, 


so  we  have 


=  ex  +8  e 
Cu  Cu  PCu 


r2X 


(144) 


T2  ^Cu  +  ^2e 


rox 


(145) 


But  from  equation  (134) 


BCu  =  @2  (1  +  a2r2)  (146) 

a2  =  0.615879 

r 2  =  -9.036047785 
a2r2  =  -5.565112074 


1 


+  a2r2  = 


-4.565112074 
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so  that 

TCu  =  “Cu  '  4.565U207432er2x 

t2  ■  “cu  +  e2er2X 

When  x  =  0 , 

T°2  =  64.3398° F 

=  66 . 1607°  F 
Cu 

T°n  -  T?n  =  1.8209  =  -5 . 565112074$ 0 
Cu  2  2 

32  =  -0.3271991607 

-4 . 5651120743 2  =  1.493700839 

=  64.3398  +  0.3272 
Cu 

ck  =  64. 6670°  F. 

Cu 

Finally , 

=  64.6670  +  1.4937  er2x 
Cu 

r2x 


(147) 

(148) 


(149) 


T  =  64.6670  -  0.32720  e 


(150) 


' 


. 


- 
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Table  12  lists  the  various  constants  needed  to  solve  the 
differential  equations.  I  have  listed  the  pure  helium  and  pure 
nitrogen  streams  for  the  low  temperature  end  of  the  interchanger 
as  the  solutions  are  of  the  same  form. 

Table  13  gives  the  equations  for  r  and  the  various  negative 
roots  resulting  therefrom. 

Table  14  gives  the  temperatures  of  the  gas  streams  and  the 
copper  tube  wall  surface  as  a  function  of  interchanger  length. 

Table  15  gives  the  additional  length  of  interchanger  tubing 
needed  at  the  warm  end  of  the  interchanger  to  reduce  the  tempera¬ 
ture  of  the  copper  tube  wall  to  within  0.01°F  of  its  final  value. 

It  is  clear  from  table  15  that  if  we  isolate  the  31  tubes  of 
the  mixture  stream  from  the  helium  and  nitrogen  streams  and  allow 
a  foot  additional  length,  we  can  be  confident  that  the  tube  wall 
surface  will  measure  the  temperature  of  the  gas  to  within  0.01° F. 

ESTIMATION  OF  ADDITIONAL  INTERCHANGER  TUBING  NEEDED  TO  HAVE  THE 

COPPER  TUBE  WALL  INDICATE  THE  FINAL  EQUILIBRIUM  TEMPERATURE 

OF  THE  PURE  HELIUM  AND  NITROGEN  STREAMS  TO  WITHIN  0.01° F 

When  streams  1  and  3  come  out  of  the  interchanger  at  the  cold 
end,  they  will  not  be  at  exactly  the  same  temperature.  Also,  the 
copper  tube  wall  surface  will  be  at  a  temperature  between  streams 
1  and  3  and  the  mixture  stream,  stream  2.  We  then  ask  the  question, 
"How  much  additional  length  of  tubing  is  required  to  have  the  tem¬ 
perature  of  the  tube  wall  within  0.01° F  of  the  final  equilibrium 


95X1 

« 

- 


* 
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TABLE 

12.  -  Constants 

for  the 

2 

equation  kAa£r 

+  kAr  -  = 

0 

ream 

Gas 

SCFM 

1°  ,°F 

a 

b 

kA 

2 

He 

20 

64.8935 

0.406337 

15.3938 

0.917517 

Cu 

66.1665 

2 

He-N0 (75-25) 

80/3 

64.7527 

0.516957 

22.8335 

0.917517 

Cu 

Z 

66.3089 

2 

He-N0 (50-50) 

40 

64.3398 

0.615879 

37.8481 

0.917517 

Cu 

Z 

66.1607 

2 

He-N0 (25-75) 

80/3 

64.3430 

0.605252 

27.7581 

0.917517 

Cu 

Z 

66.1633 

2 

N 

20 

64.1402 

0.599647 

22.8366 

0.917517 

Cu 

2 

65.9757 

1 

He 

20 

35.1065 

0.585280 

15.3938' 

0.710336 

Cu 

33.2729 

3 

N_ 

20 

35.8598 

0.661329 

22.8366 

0.710336 

Cu 

2 

33.8355 

1 

V*  ' 
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TABLE  13.  -  Equations  for  r  and  roots  for  the  various  calculations 

carried  out 


Stream 

Ga  s 

SCFM 

2 

He 

20 

2 

He-N2(75-25) 

80/ 

2 

He-N2(50-50) 

40 

2 

He-N2(25-75) 

80/ 

2 

N2 

20 

1 

He 

20 

3 

N2 

20 

r  +  2 . 4610114266r  -  41.290044434  =  0 
r2  =  -7.772996735 


r  +  1 . 934396865r  -  48.13976288  =  0 
r2  =  -7.97257035 

r2  +  1 . 623695564r  -  66.97836887  =  0 
r2  =  -9.036047785 

r2  +  1 . 652204371r  -  49.98496394  =  0 
r2  =  -7.94420648 

r2  +  1 . 6676477994r  -  41.5070301  =  0 
r2  =  -7.33015305 

r2  +  1 . 7085839255r  -  37.026983333  =  0 
r2  =  -6.998948025 

r2  +  1 . 5121066823r  -  48.612734622  =  0 
r 2  =  -7.769208925 


T  =  at-.  + 


r  0x 


Cu  1  +  a2r2  Cu 


3.  e  2 


a  +  3P  e 

Cu  Cu 


rnX 


1 

/ 
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TABLE 

Stream 

2 

Cu 

2 

Cu 

2 

Cu 

2 

Cu 

2 

Cu 

TABLE 


Stream 

2 

2 

2 

2 

2 


^•4 *  ~  Equations  for  temperature  of  the  gas  and  the  copper  tube 

wall  as  a  function  of  interchanger  length 


Gas 

SCFM 

He 

20 

He-N2(75-25) 

80/3 

He-N2(50-50) 

40 

He-N2(25-75) 

80/3 

N2 

20 

T2  =  65.2965  -  0.4030er2x 

Tr  =  65.2965  +  0.8700er2x 
Cu 

T2  =  65.1303  -  0.3776er2x 

T_  =  65.1303  +  1.1786er2x 

Cu 

T2  =  64.6670  -  0.3272er2x 
TCu  =  64  *  6670  +  1 .4937er2x 

T2  =  64.7216  -  0.3786er2x 
TCu  =  64.7216  +  1.4417er2x 

T2  =  64.5578  -  0.4176er2x 

Tr  =  64.5578  +  1.4179er2x 
Cu 


15 .  -  Additional  length  of  interchanger  needed,  at  the  warm 

end  of  interchanger ,  to  reduce  the  temperature  of 

copper  tube  wall  to  within  0.01  F  of  its  final 

value 


Gas 

SCFM 

L,  ft. 

He 

20 

0.5745 

He-N  (75-25) 

80/3 

0.5982 

He-N^ (50-50) 

40 

0.5541 

He-N2 (25-75) 

80/3 

0.6257 

N  2 

20 

0.6759 

2 


' 
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temperature  of  streams  1  and  3?" 

Again,  to  answer  this  question,  we  must  allow  for  heat  conduction 
along  the  tube  wall.  Neglecting  radial  temperature  gradients  in  the 
copper,  we  have  as  our  differential  equations,  as  before 


al  T1 


(151) 


(152) 


kAT 


// 

Cu 


Vl  +  b3T3 


(153) 


As  solutions 


to  equations  (151) , 


T  =  ot  e 
1  1 


rx 


T  =  ^oe 
Cu  2 


rx 


^3  °i3e 


rx 


(152)  ,  and  (153)  ,  we  will  try 

(154) 

(155) 

(156) 


T 


1 


cx^re 


rx 


(157) 


rx 

c*2re 


(158) 


rx 

Qf3re 


(159) 


Then 


' 


65 


T//  ^  2  rx 

T  Cu  =  V  e 

(160) 

Substituting  in  equation  (151)  ,  we  have 


\  rX  rX 
(0^2  -  c^)  e  =  ac^re 

(161) 

and  in  order  for  equation  (161)  to  be  true  for  all  x,  it  is  necessary 
and  sufficient  that 


a 2  =  (1  +  a^r) 

(162) 

Substituting  in  equation  (152)  ,  we  have 


,  N  rx  rx 

(c*2  -  e  =  a3a3re 

(163) 

and  in  order  for  equation  (163)  to  be  true  for  all  x,  it  is  necessary 


and  sufficient  that 

^2  =  a3  ^  +  a3r^ 

(164) 

Substituting  in  equation  (153)  ,  we  have 


i  a  2  rx  N  rx 

kA  ct^r  e  =  (b^c^  +  b^ayre 

(165) 

and  in  order  for  equation  (165)  to  be  true  for  all  x,  it  is  necessary 
and  sufficient  that 


2 

kAoc^r  =  (b  +  b^3)r 

(166) 

(166) 


/ 


. 
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Substituting  for  oc ^  from  equation  (162)  and  for  fr 
equation  (164)  ,  we  have 


om 


kAo^r' 


=  ( 


1  + 


+  1  +  a  r)a2r 

3 


(167) 


and  equation  (167)  will  be  true  for  all  x,  if, 


kAr  (1  +  a^r)  (1  +  a^r)  -  b^  (1  + 


a3r)  - 


b3(l 


air;  j 


=  u 


Equation  (168)  is  a  fourth  degree  equation  in  r.  There  are  four 
roots  to  this  equation.  Let  these  roots  be  r^,  r^  ,  r^,  and  r^ .  Ob¬ 
viously,  r  =  0  is  one  of  the  roots.  Let  r^  =  0.  Then  the  other  three 
roots  are  given  by  the  roots  of  the  cubic  equation: 


kAr  (1  +  a^r)  (1  +  a^r)  -  b^  (1  +  a^r)  -  b^  (1  +  a^r)  =  0  (169) 

or 

kAa^a^r2  +  kA  (a^  +  a^)r2  +  (kA  -  b^a^  -  b^a^)r  -  b^  -  b^  =  0.  (170) 


The  cross-sectional  area  available  for  heat  conduction,  along 
the  tube  for  a  3/16-inch  O.D.,  0.042-inch  wall,  per  tube,  is 


(O.D.)2  - 


4 


(I .D . ) 2 


O.D. 


3/  16" 


0.1875" 


' 

. 


I.D.  =  0.1875  -  2  X  . 042"  =  .1035" 


=  0.0087890625  in2 


(UP 


0.0026780625  in2 


IQ/J-J-.  J-  =  0.0061110000 

4 


TT  (O.D.)2  -  (I.D.)2  j 


=  0.0191982727  in2 


For  stream  1  there  will  be  five  tubes  and  for  stream  3  there 
will  be  19  tubes,  making  a  total  of  24  tubes  available  for  heat  con 
duction.  Total  available  cross-section  area  available  for  heat  con 
duct  ion  is  then 

A  =  24  x  ■J..1.9|982727_  =  0>0031997121  ft2 

144 

For  the  thermal  conductivity  of  copper,  we  have 
k  =  222  Btu  hr”1  ft”2  (deg  F)”1ft. 


Then 


kA  =  0.710336 


t 

.  ' 
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with  a 


i 


=  0.585280 


b 


1 


15.3938 


a3  =  0.661329 


b3  =  22.8366 


a  a  =  0.3870626371 


bxa3  =  10.18036636 


a  +  a3  =  1.246609 


b  a  =  13.36580525 


kAa  a3  =  0.2749445254 
kA  (a  +  a  )  =  0.8855112506 
kA  -  b^a3  "  ^3ai  =  “  22.83583561 

"bl  "  b3  =  -38-2304 
Then  equation  (170)  becomes 

0.2749445254  r3  +  0.8855112506  r2  -  22.83583561  r  -  38.2304  =  0  (171) 

We  divide  by  0.2749445254,  and  we  have 

r3  +  3.220690609r2  -  83.05615679r  -  139.0476859  =  0.  (172) 

The  three  roots  of  this  equation  are: 


r 


2 


-1.623455805 


(173) 


r 


3 


(174) 


-10.087696295 


>. 


. 
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r4  =  +8.490461495  (175) 


The  solutions  to  our  differential  equations  (151) , 

(152),  and 

(153)  are 


=  a^erlx  +  gjer2x  +  y1er3x  +  §^er4x 

(176) 

TCu  -  «2eriX  +  @2er2*  +  Y2er3*  +  ^ 

(177) 

T3  =  «3erlx  +  +  ^'3*  +  5^*4* 

(178) 

For  an  interchanger  of  infinite  length, 


00  CO  00 

T  =  T  =5  T 

1  3  Cu* 

(179) 

This  boundary  condition  indicates  that 


§!  =  52  =  §3  -  0. 

(180) 

Also  with  r  =  0 , 


*1  =  *2  =  a3 

(181) 

Equation  (162)  leads  to 


e2  =  (1  +  aLr2) 

(182) 

Y2  =  yl  (1  +  a1r3) 

1— * 

00 
00 
' — <• 

and 


' 


<  ' . 
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while  equation  (164)  leads  to 


$2  =  33  (1  +  a3r2) 


(184) 


and 


Y2  =  v3  (1  +  a3r3) 


(185) 


With  a^  =  0.585280  and  a^  =  0.661329,  we  have 


alr2  =  -°-9501762l355;  1  +  air2  =  ^ • 049782378645 


alr3  =  -5-904126890;  1  +  a1r3  =  -4.904126890 


a3r2  =  -1.07363840407;  1  +  &  r  2  =  -0.07363840407 


a3r3  =  -6.671286106;  1  +  a  r  =  -5.671286106 


Then 


3  =  20 . 0707347083 2 


33  =  -13. 5798706232 

y  =  -0.2039098952y2 

y3  =  -0.1763268474y2 


Then  our  solutions  for  equations  (176)  ,  (177) ,  and  (178)  become 


T  =  a2  +  20.07073470832er2X  -  0 . 2039098952Y2er3X  (186) 


. 


+  ' 
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TCu  =  “2  +  +  Y^3*  '  (187) 

) 

T3  =  a2  -  13.5798706232er2x  -  0 . 17632684 74y2er3x  (188) 

When  x  =  0 , 

f3.  =  35 . 8399°  F ,  ifl  =  36.0805°F,  T?n  =  33.9782°F 
1  3  Cu 


so  that 

35.8399  =  c*2  +  20 . 0707347083 2  -  0 . 2039098952y2  (189) 

33.9782  =  +  32  +  y2  (190) 

36.0805  =  c*2  -  13 . 5798706232  -  0 . 1763268474y2  (191) 

or 

1.8617  =  19 .0707347083 2  -  1 . 203909895y2  (192) 

2.1023  =  -14 . 579870623 2  -  1 . 1763268474y2  (193) 

0.09762078013=  32  -  0 . 06312865831y2  (194) 

0.1441919517  =  -3  -  0.08068156969y2  (195) 

0.2418127318  =  -0 . 1438102280y2  •  (196) 


* 
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Y2  =  -1.6814710272  (197) 

$2  =  0.097620780127  -  .1061490099216 

P2  =  -0.008528229792  (198) 

C*2  =  33.9782  +  1.681471027  +  0.008528229792 


c*2  =  35. 66820°  F 

Then , 

3  =  -0.1711678377 

33  =  +0.1158122572 
y  =  +0.3428685809 
Y3  =  +0.2964884852 

We  then  have  as  our  solutions 

T  =  35.66820  -  0.17117er2X  +  0.34287er3X 

T  =  35.66820  -  0.00853er2X  -  1.68147er3X 

Cu 

T3  =  35.66820  +  0.11581er2x  +  0.29649er3X 

with 


(199) 


(200) 

(201) 

(202) 


r2  =  -1.6234558 


. 

.  .1  •>  i  •.» 


I  r 


. 
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and 

r3  =  -10.087696 

Solving  equation  (201)  for  x,  when  T  is  within  0.01°F  of 

U. 

its  final  value,  we  obtain 

x  =  0.5506  feet. 


Then  we  have 

<Vx=0.5506  '  35 .5995°  F  (203) 

<Vx=0.5506  =  35 . 6382°  F  (204) 

(T3>x=0.5506  "  35 • 7 167°  F  (205) 

This  calculation  indicates  that  just  about  6-1/2  inches  of 
interchanger  length  are  needed  to  bring  the  copper  tube  surface  bet¬ 
ween  the  nitrogen  and  helium  streams  to  within  0.01° F  of  its  final 
value.  The  streams  at  this  point  are  about  0.12°F  apart  in  tempera¬ 
ture.  Thus,  we  do  not  have  to  worry  over  the  fact  that  streams  1  and 
3  will  be  at  different  temperatures  when  they  come  out  of  the  inter¬ 
changer  . 

In  table  16,  J  list  the  various  constants  for  additional  calcu¬ 
lations,  carried  out  in  this  section. 

In  table  17,  I  list  the  equations  for  r  and  the  roots  for  the 


calculations  carried  out  in  this  section. 
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In  table  18,  I  give  the  equations  for  the  temperatures  of  the 
gases  and  the  copper  tube  wall  as  a  function  of  interchanger  length. 

In  table  19,  I  give  the  additional  length  of  interchanger  needed, 
at  the  cold  end  of  the  interchanger,  to  reduce  the  temperature  of  the 
copper  tube  wall  to  within  0.01CF  of  its  final  value. 

Table  19  shows  that  if  we  allow  about  a  foot  of  additional  tubing 
length,  we  will  measure  a  temperature  of  the  tube  wall  which  will  be 
well  within  0.01°F  of  the  equilibrated  temperature  of  streams  1  and  3. 


■ 
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TABLE 

16.  - 

Constants 

kAa^a^r^ 

for  the  equation 

+  kA(a^  +  a^)r^  + 

(kA  -  b^a^  ■ 

■  b3aL)r  -  b]L  • 

■  b3  =  0 

Stream 

Gas 

SCFM 

T°  ,°  F 

a 

b 

kA 

1 

He 

20 

35.5505 

0.585280 

15.3938 

0.710336 

3 

N2 

20/3 

35.3765 

0.530878 

7.6122 

Cu 

33.6720 

1 

He 

20 

35.8399 

0.585280 

15.3938 

0.710336 

3 

N2 

20 

36.0805 

0.661329 

22.8366 

Cu 

33.9782 

1 

He 

20/3 

35.4078 

0.469829 

5.1313 

0.710336 

3 

N2 

20 

36.0102 

0.661329 

22.8366 

Cu 

z 

33.9241 

TABLE 

17.  - 

Equat ions 

for  r  and  roots 

for  the  various  calculations 

carried 

out 

Stream 

Gas 

SCFM 

1 

He 

20 

0.01852035843r3  +  0 .06652986652r2 

3 

N2 

20/3 

-  r  -  1.930493154  =  0 

r2  =  -1.821673673 

r3  =  -8.501313605 

1 

He 

20 

0.01204004662r3  +  0 .03877726508r2 

3 

N2 

20 

-  r  -  1.6741406207  =  0 

r2  =  -1.623455805 

r3  =  -10.087696295 

1 

He 

20/3 

0.01645559049r3  +  0 .059907246r2 

3 

N2 

20 

-  r  -  2.085221127  =  0 

r2  =  -1.978175358 
r3  =  -8.877854985 
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TABLE  18 .  -  Equations  for  temperatures  of  the  gases  and  the  copper  tube 

wall  as  a  function  of  interchanger  length 


Stream  Gas 

1  He 

Cu 

1  He 

3  N2 

Cu 


SCFM 

20  T  =  34.6589  +  0.4476erx 

Tn  =  34.6589  -  1.38598erx 
Cu 

(r  =  -6.998948025) 

20  T  =  35.11513  +  0.07360er2x  +  0.36177er3x 

20/3  T3  =  35.11513  -  0.14802er2x  +  0.40939er3x 

T  =  35.11513  -  0.00487er2x-  1.43826er3X 
Cu 


1  He  20 

3  N2  20 

Cu 


T  =  35.6682  -  0.17117er2x  +  0.34287er3x 

T3  =  35.6682  +  0.11581er2x  +  0.29649er3X 

T_  =  35.6682  -  0.00853er2x  -  1.68147er3x 
Cu 


1  He 

3  N2 

Cu 


20/3  T  =  35.54255  -  0.63108er2x  +  0.49633er3x 

20  T3  =  35.54255  +  0.14454eF2x  +  0.32311er3x 

Tr  =  35.54255  -  0.04455er2x  -  1.57390er3x 
Cu 


3  N2  20 

Cu 


r  y 

T3  =  35.4658  +  0 . 3940e 

ry 

Tn  =  35.4658  -  1.63031e 
Cu 


(r  =  -7.769208925) 
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TABLE  19.  -  Additional  length  of  interchanger,  needed  at  the  cold  end  of 
the  interchanger,  to  reduce  the  temperature  of  copper  tube 

wall  to  within  0,01°F  of  its  final  value 


Stream 

Gas 

SCFM 

L,  ft 

1 

He 

20 

0.7046 

1 

He 

20 

0.6052 

3 

N2 

20/3 

1 

He 

20 

0.5506 

3 

N2 

20 

1 

He 

20/3 

0.8157 

3 

N2 

20 

3 

N2 

20 

0.6557 

' 
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